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Introduction 

A project on droplet combustion, pursued joindy with F.L. Dryer of Princeton University, 
has now been in progress for many years. The project involves experiments on the burning of 
single droplets in various atmospheres, mainly at normal atmospheric pressure and below, 
performed in drop towers (e.g., ref. 1) and designed to be performed aboard space-based 
platforms such as the Space Shuttle or the Space Station and currently manifest for Spacelab in 
the MSL-1 flight of the Space Shuttle in April of 1997. It also involves numerical computations 
on droplet burning, performed mainly at Princeton, and asymptotic analyses of droplet burning, 
performed mainly at UCSD. The focus of the studies rests primarily on time-dependent droplet- 
burning characteristics and on extinction phenomena. 

The presentation to be given here concerns the recent research on application of 
asymptotic methods to investigation of the flame structure and extinction of alcohol droplets. 
These theoretical studies are relevant to the second of the proposed space-flight tests and are 
currently investigating the extent to which combustion of alcohols can be described by four-step 
reduced chemistry similar to that which has achieved a good degree of success for alkane flames. 
These studies have progressed to a point at which a number of definite conclusions can now be 
stated. These conclusions and the reasoning that led to them are outlined here. 

Asymptotic Descriptions of Alcohol Diffusion Flames 


A mixture fraction Z, defined to be zero in the oxidizer and unity in the fuel, serves as an 
attractive independent variable for analyses of diffusion-flame structures because its adoption 
enables the problem to be formulated in a universal manner that does not depend on the specific 
geometrical configuration. In the mixture-fraction formulation, the simplest asymptotic 
description is the Burke-Schumann reaction-sheet structure. This limit corresponds to an infinite 
Damkohler number for one-step, irreversible chemistry and provides the outer, transport-zone 
solutions to which the reaction-zone solutions are matched. In the reaction zone there is a one- 
parameter family of solutions, the parameter being Xst> the value at the stoichiometric mixture 
fraction of the scalar dissipation rate 

X = 2D| V Z | 2 , (1) 


where D is the diffusion coefficient of the outer solutions. As Xst increases, the peak temperature 
decreases, and eventually an abrupt extinction is approached, at a critical value. For droplet 
burning, it is found that 
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where di is the droplet diameter, and B is the transfer number for the droplet combustion 
process, given in an approximation of constant specific heat c p by 

B = {cJT„-n + 7QY 0 J{inW 0 $IL, (3) 

where T, W, and Y denote temperature, molecular weight, and mass fraction, the subscript °° 
identifies the ambient atmosphere, L is the energy required to vaporize a unit mass of fuel, Q 
stands for the heat released per mole of gaseous fuel in combustion, and n represents the number 
of carbon atoms in the alcohol. By use of equation (2) at Z = Z st , the droplet diameter at 
extinction die, can be obtained from Xste, the value of %st at extinction, which has been calculated 
without reference to the specific geometrical configuration. The calculation employs the reduced 
chemistry described below. 

Reduced Chemistry for Methanol 

The starting chemical-kinetic mechanism and the associated rate parameters for methanol 
are listed in Table 1. The elementary steps and rate parameters given here differ somewhat from 
those of other studies, since they are based on an evaluation of rate data performed specifically 
for the present program. In Table 1, CH 2 OH stands for the Boltzmann-equilibrated isomers 
CH 3 O and CH 2 OH, and the last two sets of rate parameters refer to this mixture. The table 
represents a minimal set, in the sense that all of the elementary steps shown there are essential. 

Introduction of steady-state approximations for all radicals except the important H atom 
results in a four-step reduced mechanism for the alcohols that can be written as 


CJH^ + jOH + yH + [p - (n + 1 + £y)] H 2 0 + i[3 n +±y 
- P (1 + a)] 0 2 -» P (H 2 + aCO) + (n - a|3) C0 2 

(I") 

CO + H 2 0 t^C0 2 + h 2 

an 

2H + M — » H 2 + M 

OH") 

3H 2 + 0 2 2H + 2H 2 0 , (IV”) 

where a is the number of CO molecules produced per H 2 molecule produced, (3 is the number of 


H 2 molecules produced per fuel molecule consumed, and y, a quantity of critical importance, is 
the number of radicals destroyed for each fuel molecule consumed. For methanol, this four-step 
description involves introducing steady-state approximations for CH 2 OH, CH 2 O, CHO, HO 2 , 
OH and O in Table 1 and results (ref. 2) in y = 2, giving a reduced mechanism that can be written 
as CH 3 OH + 2H— > 3H 2 + CO, CO + H 2 0 t;C0 2 + H 2 , 2H + M H 2 + M and 3H 2 + 0 2 t; 
2H + 2H 2 0, with corresponding rates *uC/tCh, (£ 9 /^ 3 )Ch(CcoCh 2 o/Ch 2 - Cco 2 ^ 3 /^ 9 X 
fcCAfCHCq, and *i/Ch[Q) 2 -(ChCh 0 /Cj AT^ATf)], respectively, where C/denotes the 
concentration of species *, subscript P representing the fuel, CH 3 OH, and the k's and K's are 
specific reaction-rate constants and equilibrium constants. Here approximations of partial 
equilibria have been introduced for steps 2 and 3, and the third-body concentration Cm is 
weighted with the efficiencies given in Table 1 , assuming that only N 2 , C0 2 and H 2 0 are present 
in significant concentrations in the reaction zone. 

Full numerical computations have shown that water-gas equilibrium is an excellent 
approximation in the reaction zones of these flames, much more accurate than for hydrocarbon 


HO 


flames, for example. Therefore, the second of the four steps is put into equilibrium, giving Cco 
= aCH2> where a = (Cco 2 ^j) / (Ch 2 o^s>). There results a three-step mechanism, which can be 
written as 


CH3OH + 2H + CO, -> + ocCO) + H,0 , 

a*) 

2H + «“ 1 °°» - a+T^ Hi + aC0 ) + a? 1 H *° 

(IT) 

and 
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with rates ^hQtCh, A^MCnCc^and kijCulQ^- C^C^q / (C^ATjAT^!)], respectively. 

If, in addition, a steady-state approximation for H is introduced, then there results the 
two-step reduced mechanism 

ch 3 oh ♦ 0 2 -» ( \ + + 2 “) HP ♦ ( , ] J C0 2 - ( , I J (H 2 + aCO) . 

(i) 

(iiJCa2+"a>)+(i)02-»(iiJ H J p + (i« o ) c 0 1 . 

(n) 

The steady-state approximation for H results in 


c H =(^,a:2) , ' 2 ^c^ 2 (c| i ' j 2 /c„ i o)[i -^ 14 c f )/(* 1/ c 0! )] 1 «, 

(4) 

where IcsCm Cq has been neglected in comparison with k^Cp and fci/Cc>2 because it is quite small 
in this expression at the pressures considered. The rates of steps I and n are found by reduction 
of the preceding rates to be 


- k h » %/ - 2^C m CoC h , (5) 

in which equation (4) can be used for Oh- This serves to relate the rates to the temperature and 
the concentrations of the species appearing in steps I and II, as needed in the analysis. In the 
three-step mechanism, equation (4) is discarded because Ch becomes an additional dependent 
variable described by an independent differential equation. 

The Flame Structure 

The flame structure is analyzed by rate-ratio asymptotics (RRA). In the RRA description, 
there is a thin fuel-consumption zone on the rich side of a broader but still thin layer of H2 and 
CO oxidation. This structure is illustrated in figure 1, where 8 and e are small parameters that 
measure, in the mixture-fraction coordinate, the thicknesses of the fuel-consumption zone and of 
the oxidation layer, respectively. The fuel-consumption zone, having peak temperature T°, is 
located at Z = Zo, a value displaced from Z s , by the small amount e£o, where £0 is calculated in 
the analysis. This structure is quite similar to that previously determined in the present program 
for alkane flames, such as heptane; essentially the only difference is the chemistry that occurs in 
the fuel-consumption zone. The methanol flame is in fact simpler than the alkane flames in that 
water-gas equilibrium is a much better approximation, so that a three-step description suffices. 

The RRA analysis yields the peak temperature as a function of jq 1 , as shown in figure 2. 
These results are essentially the same for both two-step and three-step descriptions. The 
minimum value of yr ( l identifies the extinction condition. From these minima, equation (2) can 


ill 


be used to calculate extinction diameters as functions of the oxygen mole fraction in the 
atmosphere, as shown in figure 3 for three different values of the atmospheric pressure p. Figure 
3 applies for pure methanol droplets that have not absorbed water during burning. In fact, water 
absorption is significant and modifies extinction conditions, as indicated below. 

Effects of Water Dissolution 

Water is produced in the flame and dissolves in the liquid fuel as it bums. We have 
completed an analysis of this effect, treating the gas phase as quasisteady. Two limiting cases 
were analyzed, namely time-dependent liquid-phase diffusion and perfect liquid-phase mixing. 
Figure 4 shows results of these analyses for the square of the droplet diameter as a function of 
time. 

It is seen from figure 4 that with time-dependent liquid-phase diffusion the d 2 plots are 
nearly straight lines. This is because the low diffusion coefficient for water in methanol leads to 
a relatively small amount of water absorption; the water is contained only in a very thin boundary 
layer at the liquid surface. On the other hand, with perfect liquid-phase mixing, the d 2 plots are 
curved, illustrating the reduced vaporization rate as water builds up inside the liquid. These 
curved lines are seen to be in better agreement with experiment. It is therefore concluded that the 
experiments involved an appreciable amount of liquid-phase mixing. This same conclusion was 
reached independently by Dryer and co-workers, on the basis of numerical integrations. The 
agreement in the conclusion from the two different approaches is a strong indication that the 
conclusion is correct. 

The water absorption also influences flame extinction, a further point of agreement with 
Dryer et al. In the present approach, since the gas phase is quasisteady, calculating extinction 
conditions entails calculating extinction diameters for droplets consisting of water-methanol 
mixtures with different percentages of water. For example, figure 5 shows predicted extinction 
diameters for droplets consisting of 50% water and 50% methanol on a mass basis. These 
diameters are seen to be appreciably larger than those in figure 3. Figure 6 illustrates how to 
determine extinction conditions with time-dependent water absorption. Basically, the droplet 
diameter and the extinction diameter are plotted as functions of time, and extinction occurs when 
these two diameters become equal. The resulting extinction predictions are in good agreement 
with experiment. 

Future Plans 

The research by asymptotics on methanol droplet combustion is essentially completed. It 
is intended to spend a little time looking at the kinetics of fuel consumption for higher alcohols, 
to determine the associated values of y ; the kinetic information needed to do a reasonable job of 
this has just become available, and so it is a timely investigation. However, most of the future 
research will involve studying heptane droplet combustion in helium-oxygen atmospheres. This 
is the system that will be employed in the first flight, MSL-1, and data on this system already are 
available from preliminary drop-tower tests. Most of the future work, for the next year or two at 
least, will therefore concentrate on this system. Data on histories of droplet and flame diameters 
and on extinction conditions will be reduced and compared with our previous theoretical 
predictions. Data on intensity profiles of emissions in OH bands also will be addressed, to 
determine what information that data can provide on flame structure, to be compared with theory. 
The general intent will be to improve understanding of heptane droplet combustion. 
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Table 1 

The Starting Chemical-Kinetic Mechanism 
and the Associated Rate Parameters 
in the form k = AT' exp (— E/R°T ) 


No.° Reactions 


A* 


E* 


If 

H + O3 — • OH -r 0 

3.52 xlO J ® 

-0.7 

17070 

lb 

OH + 0 -» H + Os 

1.15 xlO 14 

-0.324 

-175 

2f 

H 3 + 0 — OH + H 

5.06 xlO 4 

2.67 

6290 

2b 

OH + H-*HjtO 

2.22 xlO 4 

2.67 

4398 

3f 

Hj + OH — ► Hs 0 + H 

1.17 xlO* 

1.30 

3626 

3b 

HsO -i- H - H a + OH 

6.65 x 10* 

1.30 

4352 

4 f 

OH + OH ^ H3 O tO 

Jb = 5.46 x 10 n exp (0.00149T) 



4b 

HsO + 0 — » OH + OH 

k = 4.19 x 10 7 T 38 exp (0.003T) 



Hydroperoxvl Formation and Consumption 

~p 

H 4 O2 4 M — • HOa 4 M 

6.76 xlO" 

•1.4 

0 

6- 

HOs + H — ► OH + OH 

1.70 x 10 14 

0.0 

874 

7 

HOj 4 H — ♦ Hj 4 Oj 

4.28 x 10 J3 

0.0 

1411 

8 

HOs + OH — HsO + Os 

2.89 x 10 13 

0.0 

-497 


Ws 

tter-Gas Shift 



2f 

CO 4* OH — * CO3 4 H 

4.40 x 10° 

1.5 

-741 

9b 

COs -i- H - CO + OH 

4.97 x 10* 

1.5 

21446 

— 

Formyl Formation and Consumption 



.10 

CHO + H -» CO + Hs 

1.00 x 10 14 

0.0 

0 

11 C 

CHO + M -» CO + H + M 

2.85 x 10 14 

0.0 

16795 

12 

CH-jO + H — CHO + Hs 

2.26 x 10 10 

1.05 

3278 

Fuel Consumption and Formaldehyde Formation 

13 

CHsOH + H — CH-O + Hs 

3.00 x 1Q J * 

0.0 

0 

14 

CH3OH + H — CHsOH + Hs 

4.00 X 10 13 

0.0 

6092 


* Units: mol/lcm 3 , K, cal/mol. 

* Chaperon Efficiencies: Ns, Os: 1.0, CO: 1.9, CO3: 3.8, H3: 2.5, H3O: 12.0. 

* Chaperon Efficiencies: same as b except H3O: 16.3 
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Figure 3 The predicted extinction diameters^ of pure 
methanol droplets, burning in oxygen-nitrogen mixtures at 
300 K, as functions of the ambient oxygen mole fraction* 
for three different pressures. 



Figure 4 The square of the droplet diameter as a function 
of time for droplets initially 1 mm in diameter, burning in 
air tip = 1 atm in a helium-oxygen mixture with XO 22 » 
0.5 at p = 1 atm, and for a droplet initially 0.7 mm in 
diameter, burning in air at p « 0.25 atm; solid lines are for 
perfect liquid-phase mixing with the present theory, dashed 
lines are for time-dependent liquid-phase diffusion with the 
present theory, and points are data from the literature. 




Figure 5 The predicted extinction diameters of droplets 
composed of 50% methanol and 50% water on a mass 
basis, burning in oxygen-nitrogen mixtures at 300 K, as 
functions of the ambient oxygen mole fraction, for three 
different pressures. 


Figure 6 The droplet diameter, the extinction diameter, and 
the water flux fraction functions of time with time- 

dependent water absorption for perfect liquid-phase mixing 
illustrating how extinction conditions are determined. 
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